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I demonstrate that the amplitude of the high-energy scattering can be factorized 
in a convolution of the contributions due to fast and slow fields. The fast and slow 
fields interact by means of Wilson-line operators - infinite gauge factors ordered 
along the straight line. The resulting factorization formula gives a starting point 
for a new approach to the effective action for high-energy scattering. 

1 Introduction 

The starting point of almost every perturbative QCD calculation is a factor¬ 
ization formula of some sort. A classical example is the factorization of the 
structure functions of deep inelastic scattering into coefficient functions and 
parton densities.|^The form of factorization is dictated by process kinematics 
(for a review, sefJ). In case of deep inelastic scattering, there are two different 
scales of transverse momentum and it is therefore natural to factorize the am¬ 
plitude in the product of contributions of hard and soft parts coming from the 
regions of small and large transverse momenta, respectively. On the contrary, 
in the case of high-energy (Regge-type) processes, all the transverse momenta 
are of the same order of magnitude, but colliding particles strongly differ in 
rapidity. Consequently, it is natural to look for factorization in the rapidity 
space. 

The basic result of the paper is that the high-energy scattering amplitude 
can be factorized in a convolution of contributions due to “fast” and “slow” 
fields. To be precise, we choose a certain rapidity 770 to be a “rapidity divide” 
and we call fields with rj > rjQ fast and fields with rj < r]Q slow where rjo lies in 
the region between spectator rapidity and target rapidity. (The interpretation 
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of this fields as fast and slow is literally true only in the rest frame of the target 
but we will use this terminology for any frame). _ 

Our starting point is the operator expansion for high-energy scattering □ 
where the explicit integration over fast fields gives the coefficient functions for 
the Wilson-line operators representing the integrals over slow fields. For a 2=>2 
particle scattering in Regge limit s ^ (where m is a common mass scale 
for all other momenta in the problem t ^ p\ ^ ~ Pb ~ ~ we 

have: 

A{pa,Pb ^ p'a,p'b) = ( 1 ) 

X) / d'^xi...d'^XnC^^-"^ixi,...Xn){pB\TT:{Ui^{xi)...Ui^{x„)}\p'g) 

(As usual, s = {pA +Pb)‘^ and t = {pA —p'^)'^). Here Xi {i = 1, 2) are the trans¬ 
verse coordinates (orthogonal to both pA and p b ) and Ui{x) = If^^x)^U{x) 
where the Wilson-line operator U{x) is the gauge link ordered along the infi¬ 
nite straight line corresponding to the “rapidity divide” po. Both coefficient 
functions and matrix elements in Eq. (Q) depend on the po but this depen¬ 
dence is canceled in the physical amplitude just as the scale p (separating 
coefficient functions and matrix elements) disappears from the final results 
for structure functions in case of usual factorization. Typically, we have the 
factors ~ {g^lns/m^ — Vo) coming from the “fast” integral and the factors 
^ g'^Po coming from the “slow” integral so they combine in a usual log fac¬ 
tor g'^hiis/m'^. In the leading log approxim^ion these factors sum up into 
the BFKL pomerorBp (for a review see ref. □). Note, however, that unlike 
usual factorization, the expansion i) does not have the additional meaning 
of perturbative vs nonperturbative separation - both the coefficient functions 
and the matrix elements have perturbative and non-perturbative parts. This 
happens due to the fact that the coupling constant in a scattering processis is 
determined by the scale of transverse momenta. When we perform the usual 
factorization in hard {kj_ > p) and soft {k± < p) momenta, we calculate the 
coefficient functions perturbatively (because as{k± > p) is small) whereas the 
matrix elements are non-perturbative. Conversely, when we factorize the am¬ 
plitude in rapidity, both fast and slow parts have contributions coming from 
the regions of large and small k±. In this sense, coefficient functions and 
matrix elements enter the expansion (|^) on equal footing. We could have inte¬ 
grated first over slow fields (having the rapidities close to that of pb) and the 
expansion would have the form: 

Ais,t) d'^Xi...d'^XnD^^-'‘’^{xi,...Xn){pA\TT{Ui^{xi)...U^^{Xn)}\pA){‘^) 

In this case, the coefficient functions D are the results of integration over 
slow fields ant the matrix elements of the U operators contain only the large 
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rapidities rj > r]o. The symmetry between Eqs. (1) and (2) calls for a factor¬ 
ization formula which would have this symmetry between slow and fast fields 
in explicit form. 

Our goal is to demonstrate that one can combine the operator expansions 
( 0 ) and (|) in the following way: 

A{s,t) (fxi...(Pxn (3) 

(pAIC/“1 (a:i)... {Xn ) |p'^) (PBI(a:i)... {Xn)\ p'b) 

where = Tr(A“f7i) (A“ are the Gell-Mann matrices). It is possible to rewrite 
this factorization formula in a more visual form if we agree that operators U act 
only on states B and B' and introduce the notation Vi for the same operator 
as Ui only acting on the A and A! states: 

A{s,t) = (p^|(pB|exp(i/(i2a;E“(x)C/“(a:)) \v'j^\v'b) (4) 

In a sense, this formula amounts to writing the coefficient functions in Eq. 



Figure 1: Structure of the factorization formula. Dashed, solid, and wavy lines denote 
photons, quarks, and gluons, respectively. Wilson-line operators are denoted by dotted lines 
and the vector n gives the direction of the “rapidity divide” between fast and slow fields 

( 11 ) (or Eq. ( 0 )) as matrix_elements of Wilson-line operators. (Such an idea 
was first discussed in ref. B). Eq. (i illustrated in Fig.l is our main result 
and the rest of the paper is devoted to the derivation of this formula and the 
discussion of its possible applications. 
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2 Operator expansion for high-energy scattering 

Let us now briefly remind how to obtain the operator expansion (|^). For 
simplicity, consider the classical example of high-energy scattering of virtual 
photons with virtualities ~ — m^. 

A{s,t) = -i{0\T{j{pA)j{p'A)jiPB)j{p'B)}\0)- (5) 

where j (p) is the Fourier transform of electromagnetic current (x) multiplied 
by some suitable polarization e^(p). At high energies it is convenient to use 
the Sudakov decomposition: 


= OpPi + /3pP2 + P± (6) 

where p^ and P 2 are the light-like vectors close to pA and ps, respectively 
(Pa = Pi ~ P2P\/^^ Pb~P2~ PiP\/^)- We want to integrate over the fields 
with a > a where a is defined in such a way that the corresponding rapidity 
is pq. (In explicit form po = In ? where a = ^). The result of the integration 
will be givernby Green functions of the fast particles in slow “external” fieldfl 
(see also refill). Since the fast particle moves along a straight-line classical 
trajectory, the propagator is proportional to the straight-lin^ ordered gauge 
factor Ua. For example, when x+ > 0, p+ < 0 it has the fornfl: 


G{x,y) = i 


JdzS{z^) 


{¥- f) p/i 

27r^(x — z)'^ 


U{zi) 


27r^(z — uY 


(7) 


We use the notations z, = z^p^ and z* = z^P 2 which are essentially identical 
to the light-front coordinates z+ = Zt,ly/s, Z- = Ztjy/s. The Wilson-line 
operator U is defined as 


C7(x_l) = [oopi-I-xy, — oopi-I-xy] (8) 

where [x, y] is the straight-line ordered gauge link suspended between the points 
X and y. 


[x, y]“iVexp [ig du{x - yYA^{ux -f (1 - u)y)^ (9) 

The origin of Eq. ( 0 ) is more clear in the rest frame of the “A” pho¬ 
ton (see Fig.2). fields are approaching this quark at high speed. Due to the 
Lorentz contraction, these fields are squeezed in a shock wave located at z* = 0. 
Therefore, the propagator ( 0 ) of the quark in this shock-wave background is a 
product of three factors which reflect (i) free propagation from x to the shock 
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Figure 2: Quark propagator in the shock-wave background 


wave (ii) instantaneous interaction with the shock wave which is described by 
the operator U{z±_), and (hi) free propagation from the point of interaction z 
to the final destination y. 

The propagation of the quark-antiquark pair in the shock-wave background 
is described by the product of two propagators of Eq. (^) type which contain 
two Wilson-line factors U{z)U\z') where z' is the point where the antiquark 
crosses the shock wave. If we substitute this quark-rmtiquark propagator in 
the original expression for the amplitude (||) we obtaiiH: 

^d^xd^ze^P^--+^<i-^{T{j[x + z)j{z)})A ^ J ^I(p±,q±)Ti{U(p^)UHqz_ - p±)} 

where U(p±) is the Fourier transform of U(x±) and I(p±,q±) is the so-call^ 
“impact factor” which is a function of p±,p±-q±, and photon virtuality Ep. 
Thus, we have reproduced the leading term in the expansion (|^. (To recognize 
it, note that U(x±)U^y±) = Pexp|— fj dziUi{z±_)^ where the precise form 
of the path between points x± and ?/j_ does not matter since this is actually a 
formula for the gauge link in a pure gauge field Ui{zi_)). 

Note that formally we have obtained the operators U ordered along the 
light-like lines. Matrix elements of such operators contain divergent longitudi¬ 
nal integrations which reflect the fact that light-like gauge factor corresponds 
to a quarlemoving with speed of light (i.e., with infinite energy). As demon¬ 
strated ina, we may regularize this divergence by changing the slope of the 
supporting line: if we wish the longitudinal integration stop at 77 = 770 , we 
should order our gauge factors U along a line parallel to tt, = api + dp 2 ■ Then 
the coefficient functions in front of Wilson-line operators will contain loga- 
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rithms ^ g^lnl/a. For example, there are corrections of such type to the 
impact factor I{p,q) and if we sum them, the impact factor will be replaced 
by ^ ( 5 ^ In \/aY^ JC^I{p, q) where K. is the BFKL kernel. 

Factorization formula for high-energy scattering 

In order to understand how this expansion can be generated by the factoriza¬ 
tion formula of Eq. type we have to rederive the operator expansion in 
axial gauge A, = 0 with an additiraal condition = 0 (the existence 

of such a gauge was illustrated ii£3 by an explicit construction). It is impor¬ 
tant to note that with with power accuracy (up to corrections ~ a) our gauge 
condition may be replaced by = 0. In this gauge the coefficient functions 
are given by Feynman diagrams in the external field 

B,{x) = U,{x±)&{x,), (10) 

which is a gauge rotation of our shock wave (it is easy to see that the only 
nonzero component of the field strength tensor F,i{x) = Ui{x±) 6 {x^.) corre¬ 
sponds to shock wave). The Green functions in external field can be 
obtained from a generating functional with a source responsible for this exter¬ 
nal field. Normally, the source for given external field Afj, is just J^, = 
so in our case the only non-vanishing contribution is J^{B) = D’^Fi^,. However, 
we have a problem because the field which we try to create by this source does 
not decrease at infinity. To illustrate the problem, suppose that we use another 
light-like gauge A^, = 0 for a calculation of the propagators in the external field 
In this case, the only would-be nonzero contribution to the source term 
in the functional integral D^Fi^A* vanishes, and it looks like we do not need 
a source at all to generate the field B^- (This is of course wrong since B^ is 
not the classical solution). What it really means is that the source in this case 
lies entirely at the infinity. Indeed, when we are trying to make an external 
field A in the functional integral by the source we need to make a shift 
Aij, —>■ A^ A^ in the functional integral 

/ VAey.p {jS'(^) — ij d"‘xJ“(a::)^“^(x)} (11) 

after which the linear term D^^F^^A’^ cancels with our source term and 

the terms quadratic in A make the Green functions in the external field A. 
(Note that the classical action 5'(^) for our external field A = B (^^ vanishes). 
However, in order to reduce the linear term f cAxP^^'^Dij^A^ in the functional 
integral to the form J (PxD^^F^^A'^ {x) we need to make an integration by 
parts, and if the external field does not decrease there will be additional surface 
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terms at infinity. In our case we are trying to make the external field A = B 
so the linear term which need to be canceled by the source is 

I= Jdx^,d'^x±Fi,A^\^'Z°^^ (12) 

It comes entirely from the boundaries of integration. If we recall that in our 
case F,i(x) = Ui{x±)5{xA) we can finally rewrite the linear term as 

f d^x±Ui(x±){A^(—oop 2 + x±) — A^(oop 2 + a:±)} (13) 

The source term which we must add to the exponent in the functional integral 
to cancel the linear term after the shift is given by Eq. &) with the minus 
sign. Thus, Feynman diagrams in the external field ( p^ in the light-like gauge 
A* = 0 are generated by the functional integral 

JVAexp ^iS{A)+i j d^a:_LC/“*(a:_L)[^“(oop 2 +a^_L) —^“(—oop 2 +a;±)]| (14) 

In an arbitrary gauge the source term in the exponent in Eq. & can be 
rewritten in the form 


2i f d^x±TT{l7^(x±) fZZ„dv[-oop2 +x±,vp2 
F^:i{vp2 + xi_)[vp2 + xx, -oop2 + a:_L]} 


-a:_L] 


(15) 


Thus, we have found the generating functional for our Feynman diagrams 
in the external field dni)- However, it is easy to see (by inspection of the 
first rung of BFKL ladder diagram) that the longitudinal integrals over a in 
these diagrams will be unrestricted from below while we need the restriction 
a > a. Fortunately, we already faced that problem on the other side - in 
matrix elements of operators U and we have solved it by changing the slope 
of the supporting line. Similarly to the case of matrix elements, it can be 
demonstrated that if we want the logarithmical integrations over large a to 
stop at a = tr, we need to order the gauge factors in Eq. (H) along the same 
vector n = api + ap 2 , cf. Eq. (^). Therefore, the final form of the generating 
functional for the Feynman diagrams (with a > a cutoff) in the external field 

(0) is 


VAV'^exp ^iS{A,'i’) + i J d^x±U^^{x±)VZix±) 


(16) 


where 


V,(x±) = (17) 

dv[—oon + x±, vn + x±]nf^Ffj,i{vn + x±)[vn + x±, —oon -|- xj,] 
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and = Ti{X°'Vi) as usual. For completeness, we have added integration 
over quark fields so S'(-4, 4') is the full QCD action. 

Now we can assemble the different parts of the factorization formula (^. 
We have written down the generating functional integral for the diagrams 
with a > cr in the external fields with a < a and what remains now is to write 
down the integral over these “external” fields. Since this integral is completely 
independent of (|l^ ) we will use a different notation B and x for fh® a < a 
fields. We have: 

JvAV^V'^e'-^^^''^'>j{pA)j{p'A)j{PB)j{p'B) = ( 18 ) 

JVAV^V^e^s(^^^^j{pA)j{p’A) pBVxVx 

j {PB )j (p'b ) exp I f 11°-^ (a;±) (a:_L) | 

The operator Ui in an arbitrary gauge is given by the same formula (0) as op¬ 
erator Vi with the only difference that the gauge links and F,i are constructed 
from the fields This is our main result (^) in the functional integral repre¬ 
sentation. 

The functional integrals over A fields give logarithms of the type In 1 /cr 
while the integrals over slow B fields give powers of \n{as/m^). With loga¬ 
rithmic accuracy, they add up to g'^ In s/m^. However, there will be additional 
terms ~ g^ due to mismatch coming from the region of integration near the 
dividing point a ^ tr where the details of the cutoff in the matrix elements 
of the operators U and V become important. Therefore, one should expect 
the corrections of order of g^ to the effective action J dx±U^Vi. Still, the 
fact that the fast quark moves_along the straight line has nothing to do with 
perturbation theory (cf. ref. til); therefore it is natural to expect the non- 
perturbative generalization of the factorization formula ( p^ ) constructed from 
the same Wilson-line operators Ui and Vi (probably with some kind of non-local 
interactions between them). 


3 Effective action for high-energy scattering 


The factorization formula gives us a starting point for a new approach to the 
analysis of the high-energy effective action. Consider another rapidity t]q in 
the region between 770 and In m^/s. If we use the factorization formula (18) 
once more, this time dividing between the rapidities greater and smaller than 
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we get the expression for the amplitude (^) in the form: 


iA{s,t) = j{pA)j{p'A)3{PB)j{p'B 


(19) 



= JvAe"^^-^'^j{pA)j{p'A) J'DBA^^^'^j{pB)j{p'B) 



{x±)U^ {x±)+i J d^xj^W‘^^{x±)Y^{x±) 


(For brevity, we do not display the quark fields). In this formula operators 
Vi (made from A) fields are given by Eq. ©> the operators Ui are also 
given by Eq. (0) but constructed from C fields, and the operators Wi (made 
from C fields) and 1) (made from B fields) are aligned along the direction 
n' = a'Pi + a'p 2 corresponding to the rapidity p' (as usual, Via' j a' = p' where 
a' = rr? I sa')-. 

U.(xA)=fZ ^ dv[—oon + x±,vn + x±]n'^Ff^i{vn + x±)[vn + x±, —oon + a;_L]c 
Wi{x_L) = —oon' + x±, vn' + a;_L]n ^Ff^i{vn' + x±)[vn' + x±, —oon' + a;_L]c 

Yi{x±) = JZj dv[—oon' + x±,vn' + x±]n '''F^iivn' + x±)[vn' + xa, —oon' + x_i\b 

Thus, we have factorized the functional integral over “old” B fields into the 
product of two integrals over C and “new” B fields. 

Now, let us integrate over the C fields and write down the result in terms 
of an effective action. Formally, one obtains: 

iA{s,t)= JvAF^^^'>jipA)j{pZ (20) 

where S'efi for the rapidity interval between p and p' is defined as 

giSetf(Vi,Vi;^) _ <BxxV“-'{xx)U^{xx)+if <BxxW'‘'{xx)y;;-{xx) ^2^) 

This formula gives a rigorous-definition for the effective action for a given 
interval in rapidity (cf. ref. u). Next step would be to perform explicitly 
the integrations over the longitudinal momenta in the r.h.s. of Eq. and 
obtain the answer for the integration over our rapidity region (from p to p') 
in terms of two-dimensional theory in the transverse coordinate space which 
hopefully would give us the unitarization of the BFKL pomeron. At present, 
it is not known how to do this. One can obtain, however, a first few terms in 
the expansion of effective action in powers of Vi and Yi. The easiest way to 
do this is to expand gauge factors Ui and Wi in r.h.s. of Eq. (0 ) in powers 
of C fields and calculate the relevant perturbative diagrams (see Eig.2l_^or 
illustration, let us present a couple of first terms in the effective actionEjfJ: 

5eff = J d^xV<^\x)Y-ix)- (22) 
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Figure 3: Perturbative expansion of effective action. 
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■ In ^ (^Nc f (Px(PyV^^{x) \v?{x — y)‘^Yj‘j{y) + / cPxcPydPx'(Py' 

in ^ in ^ (^) ' in In + 


where we we use the notation V^j(x) = ^V°‘{x) etc. The first term (see 
Fig. 2a) looks like the corresponding term in the factorization formula (18) 
- only the directions of the supporting lines are now strongly different. The 
second termj^hown in Fig. 2b is the first-order expression for the reggeization 
of the gluoiiand the third term (see Fig. 2c) is the two-reggeon Lipatov’s 
HamiltoniaiO responsible for BFKL logarithms. 


4 Effective action and collision of two shock waves 

The functional integral (^) which defines the effective action is the usual 
QCD functional integral with two sources corresponding to the two colliding 
shock waves. Instead of calculation of perturbative diagrams (as it was done 
in previous section) one can use the semiclassical approach. This approach is 
relevant when the coupling constant is relatively small but the characteristic 
fields are large (in other words, when 1 but gVi ^ gYi ^1). In this case 

one can calculate the functional integral (^ ) by expansion around the new 
stationary point corresponding to the classical wave created by the collision of 
the shock waves. 

With leading log accuracy, we can replace the vector n by pi and the vector 
n' by p 2 - The classical equations for the wave created by the collision are: 

= 0 (23) 

= S{^X^)[^X,P2 + X±, -oop2 + a:_L]Viy*(x_L)[-oop2 + xx, ^x,p2 + a;_L] 
D>^F^^ = S{^x,)[^x^pi + XX, -oopi -I- a:_L]ViW(x_L)[-oopi -I- xx, fa:*p2 + a:_L] 
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where 




— ^ 


dv [— OOp2 + X ±, VP2 + X ±] F ^, i { vp2 + X ±)[ vp2 + x ±,— oop2 + x ±\, Y ^ 




dv[—oopi + x _ L,vpi + x _ i \ F , i{vpi + x±)[vpi + a;_L, —oopi + x_l], 


The r.h.s of the Eq. ( |^ ) is the first-order variational derivative of the source 
terms Jd?x±y°‘^{x 2 _)U°‘{x±_) and Jd‘^x±W°‘^{x±)Y°‘{x±) with respect to the 
gauge field. Also, as explained in Sect. 3, because our fields do not decrease at 
infinity there may be extra surface linear terms (cf. Eq. ©)■ The requirement 
of absence of such terms gives four additional equations 


J^.4.=-oo = S{2x./s)Y{xx), i"«L.=_oo = S{2x./s)V.{xx), (24) 

2 

= ^(-a:*)[oop2 -|-a:_L,-oop2 + a;_L]Fi(x_L)[-oop2 -l-a;_L,oop2 + x±] 
2 

^*i\x,=oo = ^(-a;.)[oopi -|-a:_L,-oopi -b a;_L]Ei(a;_L)[-oopi -|-a;_L,oopi + x±] 


The two sets and ( p^ ) define the classical field created by the collision of 
two shock waves. 

Unfortunately, it is not clear how to solve these equations. One can start 
with the trial field which is a simple superposition of the two shock waves (10) 


= 0, = e{x,)V^ + Q{x,)Y (25) 


and improve it by taking into account the interaction between the shock waves 
order by order (cf. ref. En) . The parameter of this expansion is the commutator 
g'^[Yi,Vk]. Moreover, it can be demonstrated that each extra commutator 
brings a factor In ^ and therefore this approach is a sort of leading logarithmic 
approximation. In the lowest nontrivial order one gets: 

= /*U[r.(--i), I4(.x)l - i « In (^1 - 

^ 16^ / ^^-^ x -ie + (x- z)^ + ie)[Yk{z±), V^{z±)] 

^ ~ if^ y ln(-a:|-b (a; - z)i-bie)[U/c(z_L),U''(z_L)] ( 26 ) 
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where a;| = is a longitudinal part of x^. These fields are obtained in the 

background-Feynman gauge. The corresponding expressions for field strength 
have the form: 




1 


^(1) _ 9 
27r2 


j^(i) ^ 9± 


j^(i) ^ 


87r2 

87r2 


dz± 


dz± 


—a;| + (a; — z)\ + ie 
1 


[Yk.V^ 


(27) 


—xf. + (a; — z)‘]_ + 


(a; - zY 


-i[Y,Vk]-[Yk,Y]) 


9^k[YJ,V^] , [Y,,Vk]-[Yk,V^ 


-X?. + (a; - z)\ + it 


Xt + ie 


dzj_ 


{x — zY 


-xl + {x- z)^ + ie 


9^k[YJ,VY [Y^,Vk]-[Yk,Vi] 


Xm — le 


X, + ie 


Let us now find the effective action. In the trivial order the only non-zero 
field strength components are F,® = S{^xYYi{x±) and = S{^x,)Vi{x±) 
so we get the familiar expression = f d'^XxV°‘''Y^. In the next order one 

has 


S'!!) = J d'^x _^ 2 j d^x±du 

^TrV* ([-oop2 + x±,up2 + x±]F,i{upi + x±)[up2 + a:_L, -oop2 + a:_L])^^^ -b 
TrF* ([— oopi -b a;_L, upi + x±]F^i{upi + a;_L)[Mpi -b a;_L, — oopi -b a;_L])^^^ j (28) 


We have seen above that the effective action contains In ^ (see Eq. (p2|)). 
With logarithmic accuracy the r.h.s of Eq. (28) reduces to 


S-li) = -^J d‘^xFY'^‘'Yx)F^l'>‘^{x) 

+ j d^a:_L2Tr[y*, Ei] ^[a:_L,-oop2 + a:_L]*^^^ - [x_l,- oopi-b (29) 


The first term contains the integral over d'^x = ^dx,dx^,d^x±. In order to 
separate the longitudinal divergencies from the infrared divergencies in the 
transverse space we will work in the d = 2 -b 2e transverse dimensions. It 
is convenient to perform at first the integral over a:* which is determined by 
a residue in the point a;* = 0. The integration over remaining light-cone 
variable x, factorizes then in the form dx,(x, or dx,/x,. This integral 
reflects our usual longitudinal logarithmic divergencies which arise from the 
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replacement of vectors n and n' in (^) by the light-like vectors pi and p 2 ■ In 
the momentum space this logatithmical divergency has the form f daja. It is 
clear that when a is close to cr (or cr') we can no longer approximate n by pi 
(or n' by p 2 )- Therefore, in the leading log accuracy this divergency should be 
replaced by In ^: 

noo 2 ^ r^' 2 (j 

/ dxm — = / da -> / da— = In— (30) 

Jo X, Jo a a a' 

The (first-order) gauge links in the second term in r.h.s. of Eq. ( |2^ ) have the 
logarithmic divergence of the same origin: 

[a:_L,-cx)pi [ dx^,— [ d^x±-^^^^^[Yk{z±),V'^{z±)] 

OTT J _QQ J \X ^)j_ 

[X_L,-CX)P 2 + ^ f dx,— ( d^x^—^^^-^\Yk[zi_),V‘^[ziy\ 

which also should be replaced by In ^. Performing the remaining integration 
over x± in the first term in r.h.s. of Eq. ( p^ we obtain the the first-order 
classical action in the form: 

• 5 '" = 1 ^ 1 "^/( 31 ) 

where 

Ltk = + Ytvt - Y^v^) (32) 

At d = 2 we have an infrared pole in S'Y> which must be cancelled by the 
corresponding divergency in the trajectory of the reggeized gluon. The gluon 
reggeization is not a classical effect in our approach - rather, it is a first quan¬ 
tum correction to our classical field (|2^). The relevant term can be obtained 
using the evolution equations for operators U from Ref. □. One gets: 


S,= 


(33) 

{Y-^xo.) - E“(2/i)) 


This expression coincides with the second term in r.h.s. of Eq. (P^) up to the 
terms proportional to higher commutators which we neglect here. 

Thus, the first-order expression for the effective action is the sum of Sd^\ 
S'(i), and S'r: 


5eff = J d^x^V^\x^)YJ^{x^)+ 


(34) 
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^ I<exr<eyAV.“(:c^) - - i'“‘(i'-^))} 

which coincide with (p^ ) up to the higher commutators. As usual, in the case 
of scattering of white objects the logarithmic infrared divergence ~ ^ cancels. 
For example, for the case of one-pomeron exchange the relevant term in the 
expansion of effective action is 


In fr J - V,-Y:^){x^) 


r(e) fdbn {\rh\/'nl ^ik \ \rbi\rmk 

(x-y)‘^l J Yi ^ 9 + F y 


y^'=y™)(yu) + 


167r^ 


In ^ / d^x±Vfix^)Y-^ix^) f d^y±dYjV,^(y±) 


vj’iv'j) 


(35) 


It is easy to see that the terms ^ ^ cancel if we project onto colorless state in 
t-channel (that is, replace by It is worth noting that in 

the two-gluon approximation the r.h.s. of the eq. ( ^^ gives the BFKL kernel. 

In conclusion let us mention that this semiclassical approach is suited 
for the study of the heavy-ion collisions. Indeed, for heavy-ion collisions the 
coupling constant may be relatively small due to high density (seet^). On the 
other hand, the fields produced by colliding ions are large so that the product 
gA is not small - which means that the Wilson-line gauge factors V and Y are 
of order of 1. In this case we need to know not only a couple of the first few 
terms in the expansion of the effective action, but the whole series. 
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